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A combined experimental and theoretical study of transition matrix elements of the 6p 2 Pj — > 
8s 2 Si/2 transition in atomic Cs is reported. Measurements of the polarization-dependent two- 
photon excitation spectrum associated with the transition were made in an approximately 200 
cm -1 range on the low frequency side of the 6s 2 Si/ 2 — > 6p 2 P3/2 resonance. The measurements 
depend parametrically on the relative transition matrix elements, but also are sensitive to far-off- 
resonance 6s 2 Si/2 — * np 2 Pj — > 8s 2 Si/2 transitions. In the past, this dependence has yielded a 
generalized sum rule, the value of which is dependent on sums of relative two-photon transition 
matrix elements. In the present case, best available determinations from other experiments are 
combined with theoretical matrix elements to extract the ratio of transition matrix elements for 
the 6p 2 Pj —> 8s 2 Si/2 (j = 1/2,3/2) transition. The resulting experimental value of 1.423(2) is in 
excellent agreement with the theoretical value, calculated using a relativistic all-order method, of 
1.425(2). 

PACS numbers: 32.80.-t, 32.70.Cs, 31.30. Jv,31.15.Md 



I. INTRODUCTION 

Although some of the earliest experiments in atomic 
physics were measurements of atomic lifetimes and os- 
cillator strengths Q, Q , precise determination of atomic 
transition matrix elements remains a demanding enter- 
prise @> IH HI El Some perspective on this may be gained 
by noting that, in spite of the development of a wide 
array of sophisticated experimental techniques, measure- 
ments depending directly on atomic transition matrix 
elements seldom have achieved a precision better than 
~ 0.5 % this bein g fo r the deeply studied alkali atoms 
H S II IS II !, SI3 13 IH El 13 13 111 Measure- 
ments made for non- alkali atoms, often motivated by the 
need for data in some other area of atomic physics re- 
search 

U 13 13 13 HI, have typically cited even lower 
precision. We have similarly been concerned with precise 
measurements in alkali atoms, and we have a continuing 
experimental program of precision measurement of rela- 
tive and absolute transition matrix elements in the one- 
electron atoms |3 13 13 13 13 ■ The main experimen- 
tal approaches have been based either on polarization- 
dependent Rayleigh and Raman scattering in atomic Cs, 
or on polarization-dependent two-photon spectroscopy 
applied to Na and Rb. 

In spite of the fact that many experimental E| > E l ITol 
13 13 13 and theoretical ^ 13 13 13 EJlHa ap- 
proaches have been applied to determination of atomic 
properties of atomic 133 Cs, including lifetimes or oscil- 
lator strengths and polarizabilities, only limited higher 
precision experimental data is available for many valence 



transition matrix elements in this atom. These studies 
have been motivated in part by a serious need for empir- 
ical data to extract more fundamental information from 
precise measurement of parity nonconservation (PNC) in 
this atom l33| . especially in light of the discrepancy re- 
ported in |34| between values of the weak charge Qw ex- 
tracted from high precision atomic physics experiment 
[H and the accepted Standard Model value. In order 
to clarify the situation, recent experimental efforts have 
concentrated on measurement of transition matrix ele- 
ments associated with the 6s - 6p and 6s - 7p multi- 
plet transitions 0, 0. However, the 6p - 7s and 7p - 
7s transitions, which are less experimentally accessible, 
make similarly important contributions. The availability 
of high-precision experimental data for any transitions 
between low-lying states of Cs is important for provid- 
ing additional information regarding the accuracy of the 
theoretical calculations in Cs which is crucial for the ac- 
curate analysis of Cs PNC experiment. 

In this paper we present results of our measurements 
of the spectral and polarization dependence of the 6s 
2 S 1/2 -> 6p 2 Pj -> 8s 2 S 1/2 (j = 1/2,3/2) two-color, two- 
photon transition in 133 Cs, where n = 6 is the dominant 
term. In an earlier report ,22j we described how such 
measurements could be interpreted in terms of a type 
of sum rule related to the scalar and vector transition 
probabilities. The sum rule is evaluated by fitting exper- 
imental polarization-dependent spectra to a generalized 
form containing, as fitting parameters, the relative two- 
photon transition matrix elements for the contributing 
transitions, allowing extraction of the ratio of the matrix 
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FIG. 1: Partial energy level diagram for atomic Cs, showing 
the excitation scheme used in the experiment. 

elements in the dominant term. The contributions from 
far-off-resonance transitions are small but significant for 
heavy atoms such as Cs and are evaluated theoretically. 
In the present report we describe that approach as ap- 
plied to the Cs 6s 2 Si/2 — ► 8s 2 S , 1 / 2 two- photon transi- 
tion. By combining the present measurements with pre- 
cisely determined 6s 2 51/2 ~~ ► 6p 2 Pj resonance transition 
matrix elements, the 6p 2 Pj — > 8s 2 Si/ 2 relative excited 
state matrix elements are determined. We point out that 
the same combined scheme can effectively be used to de- 
termine matrix elements for the 6s 2 S 1 / 2 — » 7p 2 Pj — » 7, 
8s 2 Si/2 transitions which, because of their importance 
for analysis of the precision parity nonconservation mea- 
surements in Cs. remain interesting cases to investigate 

m 

In the following sections, we briefly review our exper- 
imental approach, with particular attention to aspects 
new in the present study. This is followed by a de- 
scription of our experimental results and data analysis. 
We then describe the application of results of relativistic 
many-body calculations of the off-resonance terms, which 
allows the extraction of the desired two-photon transition 
matrix element ratios. Finally, the relative 6p 2 Pj — * 8s 
2 Si/2 relative transition matrix elements are determined 
and compared with theoretical calculations. 

II. EXPERIMENTAL APPROACH 

The general exp erimental approach is described in pre- 
vious reports [2l|, |2i| |2J, |25| , and so will only briefly be 
reviewed here. The basic experimental scheme [21] is 



illustrated in Figure 1, which contains a partial energy 
level diagram for the lowest few levels of atomic Cs. A 
block diagram of the experimental apparatus is shown 
in Figure 2. In the experiment, atoms in the 6s 2 S\/2 
ground level are excited by two-photon absorption to the 
8s 2 <Si/ 2 final level. In principle, all the intermediate 
np 2 Pj levels, including continuum terms, contribute to 
the total excitation probability 35] . In practice, and 
depending on the precision of the measurements, signifi- 
cant contributions are limited to the first few nearest-to- 
resonance terms. The first step of the excitation scheme 
is accomplished with an Ar + laser pumped Ti:Sapphire 
laser (laser 1) tuned in a several- hundred cm -1 energy 
range to the red of the 6s 2 S , 1 / 2 — > 6p 2 Pj transition, 
which has a hyperfme- weighted transition frequency of 
^3/2 = 11732.31 cm" 1 [HEH. With reference to Figure 
1, the detuning from one-photon resonance is defined as 
A = u)\ - CJ3/2, where the laser 1 frequency is u)±. Al- 
though detunings of magnitude greater than 200 cm -1 
were investigated, useful data was obtained only for | A | 
< 200 cm -1 . A Michelson-inteferometer type wavemeter, 
which has a precision of 10 -3 cm -1 , is used to determine 
u)±. The TkSapphire laser is passively stabilized with 
a thin-thick etalon combination and has a short term 
line width on the order of a few MHz. Long-term drifts, 
on a time scale longer than a typical data run of a few 
minutes, are dominated mainly by thermal and mechan- 
ical noise, and have a negligible influence on the experi- 
mental results. The average power is approximately 200 
mW. The laser 1 output is strongly linearly polarized, 
which is further purified by passing the beam through 
a Clan-Thompson prism polarizer. The beam is then 
passed through an electronically controlled liquid crystal 
retardation wave plate (LCR), which switches the linear 
polarization direction to one of two orthogonal linear po- 
larization directions. The resulting variable polarization 
beam is then directed to the Cs sample cell. The second 
excitation step is driven by the linearly polarized out- 
put from an external cavity diode laser (ECDL, laser 2), 
which generates an average power of 8 mW in a short- 
term (~ 1 s) bandwidth ~ 1 MHz. The laser has a fre- 
quency LO2 which nominally satisfies the two-photon reso- 
nance condition uj\ + u>2 = Wo> where luq is the frequency 
separation of hyperfme components of the 6s 2 S!/ 2 and 
8s 2 Si/ 2 levels. The degeneracy-weighted hyperfme av- 
eraged value of uiq = 24317.17 cm -1 |36J. The laser 2 
frequency is determined from the detuning according to 
UJ2 = - W3/2 - A. When A = 0, the excited state 
resonance frequency is u>2 — 12584.82 cm -1 36]. The 
commercial ECDL is piezoelectrically scannable over a 
range of approximately 15 GHz around the two-photon 
resonance. The diode laser output beam is made to be 
nearly collinear with that from the TkSapphire laser and 
the two beams are weakly focused with 0.5 m focal length 
lenses and overlapped in the central region of a sample 
cell assembly. 

The Cs vapor cell is an evacuated and sealed Pyrex 
cell having a length of about 7.5 cm and a diameter of 
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FIG. 2: Schematic diagram of the experimental apparatus, 
showing the layout of the main elements of the experiment. 
The ThSapphire Laser is Laser 1, and the Diode Laser is 
Laser 2. The Liquid Crystal Retarder controller is labeled as 
the LCR Controller. 



2.5 cm. Research grade windows are fused directly onto 
the cell body, leaving about a 2.0 cm diameter undis- 
torted central region for transmission of the laser beams. 
The cell is prepared on an oil-free vacuum system, and 
is evacuated to a base pressure of about 10 -8 mbar. A 
small amount of Cs metal is driven into the cell prior to 
removing it from the vacuum system. The cell tempera- 
ture is varied by placing it in a resistively heated oven, 
which allows heating to a typical temperature of 375 K, 
corresponding to a Cs vapor density of about 2.6 x 10 11 
cm -3 . The cell temperature is stabilized to about ±0.1 
K using a thermocouple attached to the coolest part of 
the cell; the thermocouple output is fed back to the oven 
power supply. 

Two-photon resonance signals are monitored by mea- 
surement of the 7p 2 Pj — > 6s 2 Si/2 cascade fluorescence 
at 455 nm and 459 nm. The fluorescence is collected at 
right angles to the laser beams by a short focal length 
field lens (~ 5 cm), which approximately collimates the 
fluorescence light on the cathode of a red-light sensitive 
photomultiplcr tube (PMT). The tube is protected from 
background light and the intense laser beams by a combi- 
nation of colored glass and narrow band interference fil- 
ters. Infrared transmitting colored glass filters mounted 
on the entrances to the oven housing further reduce back- 
ground light signals. The PMT output is amplified and 
the photon counting rate measured with a commercial 
100 MHz photon counter. Typical counting rates on two- 
photon resonance are ~ 10 3 s _1 . We point out that the 
operating temperature of 375 K turns out to be optimal 
for this experiment. At lower temperatures the fluores- 
cence signals are weak, while at elevated temperatures, 
and correspondingly larger Cs atom density, the vapor 
becomes optically thick to the 7p 2 Pj — > 6s 2 Si/ 2 sig- 
nal radiation. Then branching to the 7s 2 Si/2 and 6d 
2 Dj levels (see Figure 1) is enhanced, and signals in the 
observed decay channel are correspondingly reduced. 

The various instruments are globally controlled or 
monitored by a computerized data acquisition and in- 
strument control program. To illustrate the experimen- 



tal protocol, consider a typical experimental run, where 
laser 1 is set to a nominal frequency, which is measured 
by the wavemeter, which passes this value to the main 
control program. Sequentially, the photon-counting rate 
is measured alternately for collinear and perpendicular 
linear polarization directions of the two lasers. The re- 
tardance of the LCR is controlled with a specially writ- 
ten instrument driver that communicates with the LCR 
hardware via the computer parallel port. The frequency 
of the ECDL (laser 2) is then shifted by piezoelectrically 
scanning the ECDL cavity length; this is achieved by di- 
rect communication between the ECDL controller and 
the main experiment computer. The wavemeter reading 
is recorded again and the data cycle is repeated. This ex- 
perimental protocol is very effective for determining the 
shape of the excitation line, and for assessing the blend- 
ing of different hyperfine components of the excitation 
line shape. An alternate data-taking protocol was also 
employed. In this, once a proper setting of laser 2 was 
determined which minimized the influence of hyperfine 
blending (see discussion section) on the measured polar- 
ization, the laser 2 frequency was not scanned. Data was 
then accumulated by switching the laser polarization un- 
til sufficient statistics were obtained. 

For each data run, the experimental signal is deter- 
mined for two different states of relative linear polariza- 
tion of the excitation lasers. Although the absolute inten- 
sities of these components depend on many experimental 
factors, including the laser intensities, Cs density, and the 
sensitivity of the detection electronics, the intensity ratio 
is sensitive only to the relative polarization state of the 
lasers and the two-photon matrix elements. The main 
experimental observable for detailed analysis is then the 
linear polarization degree defined as 



Pl = 



Sh — Sj 
S\\ + Si 



(1) 



where S|| and Sj^ are the measured signal intensities 
when the laser beams are linearly polarized collinearly or 
perpendicularly, respectively. 

The nuclear spin of 133 Cs is I = 7/2, and so both the 
6s 2 Si/2 and 8s 2 Si/ 2 electronic states have hyperfine 
components of total angular momentum F = 3, 4. The 
hyperfine splitting in the ground 6s 2 Si/ 2 level is the 
International Frequency Standard|37j. which is approxi- 
mately 9.192 GHz, while the hyperfine splitting in the 8s 
2 Si/ 2 level is about one-tenth this value and is on the or- 
der of 0.9 GHz. Each of these splittings is larger than the 
one-photon Doppler width of several hundred MHz asso- 
ciated with the separate 6s 2 Si/ 2 — > 6p 2 Pj and 6p 2 Pj 
— * 8s 2 Sx/2 transitions, and so partial resolution of the 
hyperfine splitting is expected even when the two excita- 
tion laser beams are copropagating through the sample 
cell. This is illustrated in Figure 3, which shows the 

spectral and relative polarization dependence of the 
two-photon excitation rate at A = -80.37 cm -1 . The 
main contribution to the spectral width is due to Doppler 



4 



(3.4) 



(3,3) 



(4,4) (4,3) 



V 
• 


Copropagating beams 


i 


• • 


A = - 80.37 cm" 1 


• • 


• 


• S„ vs. Offset (GHz) 


• 


• 


o S^vs. Offset (GHz) 


• 
• 


si 




V V 



4 6 8 10 

Laser 2 Offset (GHz) 



14000 t 

12000 - 

10000 - 

I 8000 - 
o 

S 6000- 
4000 - 



Counterpropagating Beams 
A =-22.0 cm" 1 

• S„ vs. Offset 



: 



ooo 

• °° ° o • 

- nnnnnnnnnnno8°°° | %8noflnmpnnnor>opnonnrl 



50 100 150 200 250 

Laser 2 Offset (MHz) 



300 



FIG. 3: Experimental signal for frequency and polarization- 
dependent two-photon excitation at a detuning of -80.37 cm - 
from resonance. All four hyperfine transitions are shown, are 
notated by (F, F'), and are labeled by the vertical lines above 
each component. The laser 2 offset is measured from a con- 
venient starting point, denoted by the origin of the graph. 
Co-propagating beams. 



broadening, which is approximately double that for a 
one-photon transition. On the other hand, for counter- 
propagating beams, nearly complete cancellation of the 
Doppler width is expected, with a residual two-photon 
Dopplcr width ~ 10 MHz. A typical such scan is shown 
in Figure 4, which corresponds to the F = 4 — > F' = 4 
transition with A = -22.0 cm -1 . There it is seen that 
the spectral width of the line shape is approximately 40 
MHz, which corresponds well to the combined residual 
Dopplcr width and the spectral width of the Ti:Sapphire 
laser. It is also seen that the excitation spectrum depends 
significantly on the relative polarization state of the two 
laser beams, giving in this case a large linear polarization 
degree of approximately 0.78. 

In spite of the evidently good signal to noise ratio for 
the data in Figure 4, it proved difficult to achieve the 
desired reproducibility in the polarization measurements 
with the passively stabilized TkSapphire laser. The main 
reason for this was that measurement of the two different 
polarization states were made sequentially, and fluctua- 
tions in the laser frequency on the switching time scale 
introduced unwanted noise in the extracted polarization 
values. The reason for this is that when the lasers are 
counterpropagating through the interaction region of the 
cell, the Doppler width of the two-photon transition is 
greatly reduced, and is on the order of the natural width 
of the final 8s level. In this case, in order to obtain reli- 
able polarization measurements, the combined frequency 
drift of the two lasers would need to be smaller than a 
fraction of one MHz/s. Instead, we recorded the excita- 
tion spectrum for copropagating laser beams, for which 
the Doppler width is large (several hundred MHz), and 
for which the characteristic laser frequency drift rates of 
a few MHz/s are entirely acceptable. Note that variation 
of the polarization in the few MHz range is totally neg- 



FIG. 4: Experimental signal for frequency and polarization- 
dependent two-photon excitation at a detuning of -22.00 cm -1 
from resonance. The F = 4 — > F' = 4 transition is shown. The 
laser 2 offset is measured from a convenient starting point, 
denoted by the origin of the graph. Counterpropagating laser 
beams. 
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FIG. 5: Experimental signal for frequency and polarization- 
dependent two-photon excitation at a detuning of -167.9 cm" 
from resonance. The F = 3^F' = 4andF = 3->F' = 3 
transitions are shown. The laser 2 offset is measured from a 
convenient starting point, denoted by the origin of the graph. 
Copropagating laser beams. 



ligible. The small penalty to be paid for this is that the 
polarization values need to be determined in the wings 
of the partially blended lines. However, these determina- 
tions were made sufficiently far into the wings that the 
other hyperfine transition made nearly negligible contri- 
bution to the polarization value. In addition, the signals 
were significantly less noisy than with counterpropagat- 
ing beams, making possible consistent and repeatable 
measurements of the polarization at each detuning. A 
typical higher resolution data run is presented in Figure 
5, which corresponds to the F = 4 — > F' = 3, 4 transitions 
at A = -167.9 cm" 1 . 

Finally, we have modeled the residual effect of blending 
of the hyperfine lines on the measured polarization val- 
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ues. In the model, we approximate each hyperfine transi- 
tion line shape to be Gaussian and having the measured 
width determined by the Doppler broadening of the two- 
photon transition. Using the known relative hyperfme 
transition probabilities gives a polarization-dependent 
correction that ranges from 0.0005 to 0.0013 (polariza- 
tion values reported here can range from -1 to +1). Even 
though these corrections are well within the statistical er- 
ror associated with each measured polarization, the cor- 
rections are systematic and so are made directly to the 
measured counting rates, prior to extracting values 
from the data through Eq. (1). 



A. Analysis and Experimental Results 

The variation of the linear polarization degree with de- 
tuning A of laser 1, and for each of the four electric-dipole 
allowed hyperfme transitions, may be readily calculated 
in terms of transition matrix elements (the general ex- 
pression may be found in |35jh The dominant interme- 
diate levels for the two-photon transitions are the 6p 2 Pj 
levels, with j = 1/2, 3/2, but other np 2 Pj transitions, 
including the p-continuum, also contribute to the total 
transition probability, and play an important role in the 
experiments reported here. For the transitions when AF 
= ±1, absorption of two photons with collinear polariza- 
tion directions is forbidden, and so the linear polarization 
degree is -1, independent of detuning. However, for the 
other transition pair, where AF = 0, both fine struc- 
ture multiplet components contribute, leading to strong 
spectral variations in the linear polarization degree. The- 
oretical expressions for the intensities arc given by 
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FIG. 6: Model dependence of the variation of the linear po- 
larization degree with detuning from resonance for all four 
permitted hyperfme transitions. The matrix element ratio is 
taken to be R = 2, and contributions from far-off-resonance 
transitions are neglected. 



for the F = 4 — > F' = 4 transition. In these expres- 
sions, the overall intensity for each hyperfine transition 
is proportional to 1 ff', while R, P and Q are parame- 
ters that are pro por tional to ratios of reduced transition 
matrix elements |3!|. In the present case, 



R 



< 8s || d || 6 P3/2 >< 6 P3/2 || d || 6s > 
<8s || d || 6pi /2 >< 6pi/2 IN II 6s > : 



(6) 



where d is the electric-dipole operator. The quantities 
P and Q are given by 



I\\ = 4/o33 
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w l _ ^3/2 U\ - U>x/2 
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+ P . 
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I\\ = Io33 
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^2 — ^3/2 w 2 — Wi/2 

for the F = 3 — > F' = 3 transition and 
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Q 



(2) 



(3) 



(4) 



(5) 



P = 



E 



M, 



njk 



E 



(_l)i+Va M , 



njk 



j + 1/2 UJ k - U) n 



(7) 



(8) 



n>6,jk ■>•-'- - Jn Vi 

where the matrix element ratio is given by 



M?njk — 



< 8s || d || npj >< npj \\ d || 6s > 
<8s || d || 6 Pl/2 >< 6 Pl/2 ||d||6s>' 



(9) 



The total angular momentum j = 1/2,3/2, while k = 
1,2 labels the frequency of the absorbed photons. 

An illustration of the variation of the polarization with 
detuning A is shown in Figure 6 for each of the four hy- 
perfine transitions. In these plots we have taken R = 
2, corresponding to no relativistic modification of the re- 
duced transition dipole matrix elements, and P = Q = 0, 
which applies to the case when only the n = 6 intermedi- 
ate levels are considered. Among the critical features of 
the plots are the resonance values of the linear polariza- 
tion degree for each hyperfine transition. In the present 
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FIG. 7: Measured linear polarization degree of the two- 
photon excitation rate as a function of detuning for the F 
= 3 -> F' = 3 and F = 4 -> F' = 4 transitions. The solid 
curves represent the model curves for a matrix element ratio 
R = 2, and neglecting far-off-resonance transitions (P=Q=0). 



case, these values are V L = 0.882 for the F = 4 — > F' = 4 
transition, and P L = 0.811 for the F = 3 -> F' = 3 tran- 
sition. In addition, there is a detuning from resonance 
where the polarization is Pl = - 1.0, corresponding to 
the cases when the intensity S|| =0 for each hypcrfinc 
transition. The location of this point generally depends 
on the four main terms in Eq. (2) - (8), and on the values 
for P and Q. 

Polarization measurements have been made in an ap- 
proximately 200 cm -1 range of detunings to the low 
frequency side of the atomic 6s 2 Si/ 2 — > 6p 2 Pj res- 
onance transition. These measurements, which have a 
typical one standard-deviation uncertainty of 0.002, are 
presented in Figure 7. The solid curves in this figure are 
expanded versions of those in Figure 6. Note that the 
error bars on the experimental data points are negligi- 
ble on the scale of the figure. It can be seen from the 
curves that there is significant discrepancy between the 
measurements and the theoretical curves (with P = 0, 
Q = 0, and R = 2) for the F = 4 -> F' = 4 transition, 
for which we have taken the most extensive data. At 
the largest detuning, this amounts to about 20 standard 
deviations. The departure has two main contributions, 
these being the fact that nonresonant transitions with n 
> 6 make a significant contribution through nonzero val- 
ues of P and Q, and because the value of R generally de- 
parts from the nonrelativistic value of R — 2. In previous 
work, the three parameters P, Q and R were considered 
as constants, which is quite well justified by their weak 
dependence on detuning. Fits to the detuning-dependent 
polarization data then yielded linear relationships among 
the parameters; estimates of P and Q further allowed 
extraction of a quite precise value for the constant ma- 
trix element ratio R. Although this is found to be in- 
sufficient for the precision of the data presented here, 
we report here our values, making the same assumptions 



as in earlier work. That approach yields R=2. 1068(36) 
+ 443.2(1.2)P - 362(362)Q, where the uncertainty (in 
parenthesis) in the numerical values represents one stan- 
dard deviation. Note that although the coefficient of Q 
has an uncertainty on the order of its value, its value 
is correlated with that of R and P, and so it cannot be 
neglected. The final uncertainty in R derived from calcu- 
lated values of P and Q shows that the correction due to 
Q is on the order of the uncertainty in R. This equation 
represents the equivalent sum rule to those presented in 
earlier reports on our measurements in Rb. 

However, P and Q depend weakly on detuning, and so 
to obtain the highest precision, it is desirable to calculate 
P and Q as a function of detuning directly from the most 
accurate experimental and theoretical matrix elements. 
We use a combination of the most reliable experimental 
and theoretical values for dominant contributions with 
n < 9. For n = 6, measured resonance line (n = 6) 
matrix elements of Rafac, et al. [l^ are employed, while 
for the second resonance doublet (n = 7), we use the 
recently reported precision measurements of Vasilyev et 
al. [13j . Other values having n < 9 are obtained from 
the relativistic all-order calculations of Safronova et al. 
pjof: these are listed in Table 1. For multiplets with n 
> 9, which make only a small overall contribution to R, 
we use Dirac-Hartree-Fock matrix elements 28]; these 
values are also summarized in Table 1, along with the 
energies associated with each intermediate np level. For 
reference, it is found that P and Q, calculated with these 
values, are fit very well by second order polynomials in 
detuning A = u>\ - uj\q, and are given by P = -2.437 x 
10~ 4 + 3.257 x 10~ 9 A - 3.8847 x 10~ 12 A 2 and by Q = 
3.260 x 10~ 6 - 7.683 x 10~ 9 A + 1.645 x 10~ 13 A 2 . The 
detuning A is in vacuum cm" 1 units. We estimate the 
total uncertainty in P and Q to be on the order of 1%. 

To obtain the matrix clement ratio defined in Eq. (6), a 
nonlinear least squares fit of the polarization data of Fig. 
7 is made to the theoretical expressions, using the calcu- 
lated values of P and Q for each detuning. The result of 
the fit is R = 2.0024(24), where the dominant error in 
R comes from the statistical uncertainty in the measured 
polarization. The quality of the fit is illustrated in Fig. 8. 
In the figure, measurements are compared to the devia- 
tion of the polarization from the values obtained when R 
= 2, P = Q = 0; these are represented by the horizontal 
line passing through zero deviation. The experimental 
measurements, represented by the data points, are seen 
to be in excellent agreement with the solid curve, this 
being calculated from the fitting parameters. 

The ratio of excited-state dipole matrix elements 
R 8s - 6p = <8s||d||6p3/2>/<8s||d||6pl/2> may be ob- 
tained from Eq. (6) by combining the measured value 
of R with the high precision 6p resonance line matrix 
element ratio of 1.4074(3) measured by Rafac et al. 
[Tl|. This gives the excited state matrix element ratio 
of R8 S -6p = 1-423(2), which has the uncertainty of about 
0.15%, with roughly equal contributions from statistical 
uncertainty in the measurements reported here and from 
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FIG. 8: Deviation of the measured linear polarization degree 
of the two-photon excitation rate as a function of detuning 
for the F = 4 — > F' = 4 transition. The solid curve through 
the data points represents the best fit to the experimental 
data, while the horizontal line is the model curve result. Note 
that the vertical scale is magnified 20 times in comparison to 
Figures 6-7. 



the resonance line matrix element ratio of |TT| • Finally, 
we point out that the analysis may be reversed, if it is 
assumed that the bare transition matrix element ratio is 
calculated precisely. Then, by recognizing that the off- 
resonance terms represented by P and Q contribute an 
average of about -0.1091(37) to the value of R, the mea- 
surements determine the much smaller combination of 
dipole matrix elements in Eqs. (7) - (8) to about 3 %. 
As these terms are dominated by transitions through the 
6s - 7p and 6s - 8p multiplets, the measurements serve as 
a consistency check on existing measurements 13] at that 
level. 



III. RELATIVISTIC ALL-ORDER 
CALCULATION OF TRANSITION MATRIX 
ELEMENTS 

In order to compare the experimental result R^ s -6p — 
1.423(2) with the high-precision theoretical value we 
carry out relativistic all-order calculation of electric- 
dipole matrix elements in Cs. In particular, we calculate 
8s — np, n = 6, 7, 8, 9 electric-dipole reduced matrix 
elements using a relativistic all-order method including 
single and double (SD) excitations |28j. The resulting 
values of the 8s — 7p, 8s — 8p, and 8s — 9p electric-dipole 
matrix elements were used to determine P and Q. 

In the SD all-order method, the wave function of the 
valence electron v is represented as an expansion 



1*,) = 



Pninab&m®^n®'bQ'a 



ma mnab 
Em/n Prnv a ln a v H~~ "Ylmna PmnvaQ> m a 'liQ>aQ'V \&v )(io) 



In this equation |4>„) is the lowest-order atomic state 
wave function, which is taken to be the frozen-core Dirac- 
Hartree-Fock (DHF) wave function of a state v and a\ 
and a, are creation and annihilation operators, respec- 
tively. The indices a, b designate core electrons and in- 
dices m, n designate any states above the core. The 
equations for excitations coefficients p m a, Pmv, Pmnab, 
and Pmnva are obtained by substituting the wave func- 
tion \t„ into the many-body Schrodinger equation 



-El*,; 



(ID 



where H is the relativistic no-pair Hamiltonian 29] . The 
equations are solved iteratively until the corresponding 
correlation energy for the state v is sufficiently converged. 
The resulting excitation coefficients are then used to cal- 
culate matrix elements. The one-body matrix element of 
the operator Z given by 



(#u,|£|*„) 



is expressed in terms of excitation coefficients as 
= z wv + + ■■■ + ZW 

y/(l+N v )(l + N w ) ' 



(12) 



(13) 



where z wv is the lowest-order DHF matrix element, the 
terms Z k , k = a...t are linear or quadratic functions of 
the excitation coefficients, and normalization terms N v 
are quadratic functions of the excitation coefficients. As 
a result, certain sets of many-body perturbation theory 
(MBPT) terms are summed to all orders. This method 
is shown to yield high-accuracy results for the primary 
transition electric-dipole matrix elements in alkali-metal 
atoms [27l|23II3. The results for the reduced electric- 
dipole matrix elements for 8s — 6p, 8s — 7p, 8s — 8p, and 
8s — 9p transitions are listed in Table [I] 

Next, we investigate the effect of the correlation to the 
8s — 6p matrix elements to evaluate the uncertainty of 
their ratio. We find that the total correlation correction 
to the 8s — 6p matrix elements is small, 3%. However, it 
actually results from severe cancellations of large contri- 
butions illustrated in Table ITU where we give the break- 
down of the all-order SD calculation for both 8s — 6p 
transitions. The lowest-order (DHF) value is listed in 
the first row of Table HJ Three larger terms, Z<- a \ Z^, 
and Z^ 



Z (a) - Vz o +Vz 0* 

/ J /C, amf J 'W7nva \ / J ^raa^vr, 
ma ma 

Z {c) - Vz o +Vz o* 

m m 

% ^ — ^ 1 z mnPmwPnv) 
inn 



where 



(14) 



(15) 
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TABLE I: Energy and matrix elements (ME) used in calcu- 
lating the contributions of far-off-resonance transitions to the 
measured linear polarization spectrum. ME for n > 10 are 
obtained from Dirac-Hartree-Fock calculations. 



Relay Level Energy (cm ) 6s-npj (a.u.) 8s-npj (a.u.) 



6pl/2 


11178.2 


4.489 a 


-1.027 d 


6p3/2 


11732.4 


6.324 a 


-1.462 d 


7pl/2 


21732.4 


0.276 6 


-9.251 d 


7p3/2 


21765.7 


0.586 6 


-14.00 d 


8pl/2 


25709.1 


0.081 c 


17.71 d 


8p3/2 


25791.8 


0.218 c 


24.46 d 


9pl/2 


27637.3 


0.043 6 


1.743 d 


9p3/2 


27682.0 


0.127 6 


2.969 d 


10pl/2 


28727.1 


0.047 


0.634 


10p3/2 


28753.9 


0.114 


1.158 


llpl/2 


29403.7 


0.034 


0.348 


llp3/2 


29421.1 


0.085 


0.667 


12pl/2 


29852.9 


0.026 


0.228 


12p3/2 


29864.7 


0.067 


0.451 


13pl/2 


30166.0 


0.021 


0.165 


13p3/2 


30174.5 


0.055 


0.334 


14pl/2 


30393.2 


0.017 


0.127 


14p3/2 


30399.5 


0.046 


0.262 


15pl/2 


30563.3 


0.015 


0.102 


15p3/2 


30568.0 


0.039 


0.213 



a Rafac, et al 12] 
b Vasilyev, et al 1131 
c Safronova, et al |30| 
d This work (all order) 



TABLE II: Contributions to SD all-order electric-dipole 8s — 
6p matrix elements in Cs. Matrix elements are given in a.u. 



Contribution 


8s - 6pi/ 2 


8s - 6p 3/2 


DHF 


1.0584 


1.5145 


z (a) 


0.0173 


0.0177 


Z (c) 


-0.1122 


-0.1606 


Z (d) 


0.0938 


0.1317 


■Mother 


-0.0084 


-0.0108 


^norm 


-0.0228 


-0.0308 


Total 


1.0260 


1.4618 



are listed separately and all other terms are summed 
together as Z ot her- The total normalization correction 
■Znorm defined as 

Znorm = Z wv - [z wv + +■■■+ Z&] (16) 

is given separately. We find that two dominant terms, 
Z^ and Z( d \ nearly cancel each other and their sum 
almost exactly cancels out the term Z^ a \ 

We note that while terms Z^ and Z^ contain third- 



order terms as well as higher-order terms (see |27| for 
details) the Z^ term, being quadratic in valence single- 
excitation coefficients p mv contains only terms starting 
from fifth order. Out of the remaining terms, the largest 
contribution comes from the normalization correction. 
For these transitions, the dominant contribution to 
in the denominator of Eq. (|13fl comes from the term 

^PmvPmv (17) 
m 

This term contributes 97% to the iVg s value and 64% to 
the Nq p value. Again, this term is quadratic in single 
excitation coefficients p mv and can contain only MBPT 
terms starting from fifth order. Thus, 8s — 6p ma- 
trix elements present an interesting case with large, but 
canceling, contributions from high orders in many-body 
perturbation theory. This cancellation occurs for both 
8s — 6pi/ 2 and 8s — 6p 3 / 2 matrix elements. 

To evaluate the uncertainty in the theoretical ratio 
Rss-ep we also calculate 8s — 6^3/2 and 8s — 6pi/2 
electric-dipole matrix elements in different approxima- 
tions and estimate some omitted contributions. The 
results are summarized in Table IIIII The lowest-order 
DHF values are listed in column labeled DHF. The re- 
sults of the third-order many-body perturbation-theory 
calculation, which includes higher-order random-phase- 
approximation terms as described in M 

, are listed in col- 
umn "III" . The third-order values, which include an esti- 
mate of the omitted fourth- and higher-order Brueckner- 
orbital (see classification and formulas in |2jJ) correc- 
tions obtained by the scaling procedure described in |3l| , 
are listed in column "III SC " ■ Single-double all-order data 
from Tabled are listed in column "SD" . The results ob- 
tained by including partial contribution of the triple exci- 
tations are listed in column labeled "SDpT" . These data 
are obtained by adding a triple-excitation valence term 
to Eq. HI U| l and making corresponding corrections to cor- 
relation energy and single-valence excitation coefficient 
equations p mv as described in Refs. |27], |29|, |30j . The p mv 
excitation coefficients give rise to the largest terms, Z^ 
and Z^ d \ and the dominant part of N v . 

All-order SD and SDpT calculations include a com- 
plete third-order contribution (see j^J for a detailed 
comparison) but omit some classes of higher-order terms 
starting from fourth order. We have estimated some 
omitted correlation corrections resulting from triple and 
higher excitations using the scaling described, for exam- 
ple, in Ref. [2|j. Briefly, single-particle excitation coeffi- 
cients p mv are multiplied by the ratio of the experimental 
and theoretical correlation energies. The modified exci- 
tation coefficients are then used to re-calculate matrix el- 
ements. Such scaling estimates only certain classes of the 
omitted contributions (mainly Brueckner-orbital contri- 
butions, see Ref. for classification of the perturbation 
theory terms). The corresponding SD and SDpT scaled 
values are listed in columns "SD SC " and "SDpT sc ". Ab 
initio and scaled values are grouped together. 

The third-order matrix elements do not differ signif- 
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TABLE III: Electric-dipole 8s-6p-, reduced matrix elements in 
Cs calculated using different approximations: Dirac-Hartree- 
Fock (DHF), third-order many-body perturbation theory 
(III), single-double all-order method (SD), single-double all- 
order method including partial triple contributions (SDpT) 
and the corresponding scaled values. Rs s -a p is the ratio of 
the 8s — 6p 3 / 2 and 8s — 6pi/ 2 matrix elements. Absolute values 
of the matrix elements in a.u. are given. 



Level 




Ab initio 




Scaled 






DHF 


III 


SD 


SDpT 


IILc 


SD SC 


SDpT sc 


j=l/2 


1.0584 


1.0231 


1.0260 


1.0321 


1.0315 


1.0223 


1.0327 


j=3/2 


1.5145 


1.4519 


1.4618 


1.4709 


1.4712 


1.4556 


1.4705 


R8s-6p 


1.4309 


1.4191 


1.4247 


1.4252 


1.4262 


1.4238 


1.4240 



icantly from the all-order results indicating very accu- 
rate cancellation between large higher-order terms. The 
Brueckner- orbital correction is dominant in a third-order 
calculation for both transitions and is relatively more im- 
portant for the 8s — 6P3/2 transition, thus giving a ratio 
differing from more accurate SD all-order value by 0.6%. 
We note that the ratio of the scaled third-order values, 
which includes an estimate of higher-order Brueckner- 
orbital terms, agrees with the all-order result. As one 
can see from Table IIIII the total correlation correction 
contribution to the ratio Rg s -6p is very small, 0.4%, ow- 
ing to the cancellation of dominant terms. The difference 
between all-order data in different approximations is only 
0.1%. As the largest contributions to these matrix ele- 
ments contain valence single-excitation coefficients p mv 
the inclusion of the partial triple contributions and the 
scaling described above should give a good estimate of 
the omitted higher-order terms as both these methods 
are aimed at correcting p mv . We recommend the SD 
ab initio value Rs s -6p — 1-425(2) as the final theoreti- 
cal value for the ratio. The uncertainty is obtained by 
combining the uncertainty of the dominant terms, deter- 
mined to be 0.1% based on the spread of SD, SDpT and 
scaled values, and the total uncertainty in all other, much 
smaller, contributions taken to be 0.1% not to exceed the 
uncertainty of the dominant terms. 

IV. DISCUSSION OF RESULTS AND 
CONCLUSIONS 

Comparison of the theoretical matrix element ratio of 
R = 1.425(2), and the corresponding experimental re- 



sult of R = 1.423(2) shows excellent agreement. This 
is remarkable, given the various contributing theoreti- 
cal and experimental factors that significantly affect the 
final result in each case. It is interesting to reiterate 
that many of the two-photon transition matrix elements 
in heavy alkali atoms show significant relativistic mod- 
ification. For example, the 5s 2 Si/2 — > 5p 2 Pj — > 5d 
2 D$/2 transition in Rb shows such variations at a level of 
nearly 7 %. At the same time, other ratios, such as that 
measured in the present case, show significantly smaller 
relativistic modification, and yet in each case such modi- 
fications are generally well-described by the highest-level 
atomic structure calculations. Nevertheless, there is a 
significant sensitivity in the experimental measurements 
to other dipole transitions, and these transitions depend 
significantly on relativistic contributions. Naturally, it 
would be of interest to measure the absolute oscillator 
strengths of the individual multiplet transitions, in order 
to make direct comparison of the matrix element values 
reported here. Such measurements remain a significant 
challenge to experimental technique. Finally the excited 
state transitions in atomic Cs associated with the 7p dou- 
blet are of particular interest, as those second resonance 
line matrix elements show nearly a factor of 2 departure 
from the expected line strength ratio of 2. 

In conclusion, we have used precision two-photon po- 
larization spectroscopy to make measurements of the 
transition matrix element ratio associated with the 6s 
2 Si/2 — > np 2 Pj —* 8s 2 Si/2 transition in atomic Cs. 
The measurements are combined with other experimen- 
tal data and calculations in order to extract the ratio of 
excited-state matrix elements. The experimental value 
is found to be in excellent agreement with the ratio of 
reduced matrix elements calculated using a relativistic 
all-order method. 



Acknowledgments 



The financial support of the National Science Founda- 
tion (NSF-PHY-0099587) and Central Michigan Univer- 
sity is greatly appreciated. 



[1] F. Molish and B.P. Oehry, Radiation Trapping in Atomic 
Vapours, Clarendon Press (Oxford, 1998); Alan Corney, 
Atomic and Laser Spectroscopy, Clarendon Press (Ox- 
ford, 1998). 

[2] W.L. Wiese, M.W. Smith, and B.M. Glennon, Atomic 



Transition Probabilities, v. 1-2, Natl. Bur. Stand. (U.S.) 

Circ. No. 4 (U.S. GPO, Washington, D.C., 1966). 
[3] W.I McAlexander, E.R.I. Abraham, and R.G. Hulet, 

Phys. Rev. A 54, R5 (1996). 
[4] W.I. McAlexander, E.R.I.Abraham, N.W.M. Ritchie, 



10 



C.J. Williams, H.T.C. Stoof, and R.G. Hulet, Phys. Rev. 
A 51, R871 (1995). 

C.W. Oates, K.R. Vogel, J.L.Hall, Phys. Rev. Lett. 76, 
2866 (1996). 

U.Volz, M. Jajerus, H. Liebel, A. Schmitt, H. 

Schmoranzer, Phys. Rev. Lett, 76, 2862 (1996). 

U. Volz and H. Schmoranzer, Physica Scripta T 65, 48 

(1996). 

C.E. Tanner, A.E. Livingston, R.J. Rafac, F.G. Serpa, 
K.W. Kukla, H.G. Berry, L. Young, and C.A. Kurtz, 
Phys. Rev. Lett. 69, 2765 (1992). 

R.J. Rafac, C.E. Tanner, A.E. Livingston, K.W. Kukla, 
J.G. Berry, and C.A. Kurtz, Phys. Rev. A 50, R1976 
(1994). 

L. Young, W.T. Hill III, S.J. Sibener, S.D. Price, C.E. 
Tanner, C.E. Wieman, and S.R. Leone, Phys. Rev. A 50 
(1994). 

R.J. Rafac and C.E. Tanner, Phys. Rev. A 58, 1087 
(1998). 

R. J. Rafac, C. E. Tanner, A. E. Livingston, and H.G. 
Berry, Phys. Rev. A 60, 3648 (1999). 
A. A. Vasilyev, I.M. Savukov, M.S. Safronova, and H.G. 
Berry, Phys. Rev. A 66, 020101 (2002). 

C. R. Ekstrom, J. Schmiedmayer, M.S. Chapman, T.D. 
Hammond, D.E. Pritchard, Phys. Rev. A 51, 3883 (1995). 
K.M. Jones, PS. Julienne, PD. Lett, W.D. Phillips, 
E. Tiesinga, and C.J. Williams, Europhys, Lett. 35, 85 
(1996). 

H. Wang, J. Li, XT. Wang, C.J. Williams, PL. Gould, 
W.C. Stwalley, Phys. Rev. A 55, R1569 (1997). 
A. A. Radzig, B.M. Smirnov, Reference Data on Atoms, 
Molecules and Ions (Springer, Berlin, 1985). 
Cohen- Tannoudji, J. Dupont-Roc, G. Grynberg, Atom- 
Photon Interactions: Basic Processes and Applications 
(Wiley, New York, 1992). 

Harold J. Metcalf and Peter van der Straten, Laser Cool- 
ing and Trapping (Springer- Verlag New York, 1999). 

D. Mihalas, Stellar Atmospheres (Freeman, San Fran- 
cisco, 1970). 

A.I. Beger, M.D. Havey, R.P. Meyer, Phys. Rev. A 55, 



3780 (1997). 

[22] R.P. Meyer, A.I. Beger, and M.D. Havey, Phys. Rev. A 

55, 38780 (1997). 
[23] M.D. Havey, Phys. Lett. A 240, 219 (1998). 
[24] S.B. Bayram, M. Havey, M. Rosu, A. Sieradzan, A. Dere- 

vianko, and W.R. Johnson, Phys. Rev. A 61, 050502(R) 

(2000). 

[25] A. Markhotok, S.B. Bayram, A. Sieradzan and M.D. 

Havey, J. App. Phys. 92, 1613 (2002). 
[26] S.C. Bennett, J.L. Roberts, and C.E. Wieman, Phys. 

Rev. A 59, R16 (1999). 
[27] M.S. Safronova, Ph.D. Thesis, University of Notre Dame 

(2002). 

[28] S.A. Blunder!, W.R. Johnson, Z.W. Liu, and J. 

Sapirstcin, Phys. Rev. A 40, 2233 (1989). 
[29] S.A. Blundell, W.R. Johnson and J. Sapirstein, Phys. 

Rev. A 43, 3407 (1991). 
[30] M.S. Safronova, W.R. Johnson, and A. Derevianko, Phys. 

Rev. A 60, 4476 (1999). 
[31] W.R. Johnson and Z.W. Lui and J. Sapirstein, At. Data 

Nucl. Data Tables 40, 2233 (1989); W.R. Johnson, Z.W. 

Liu, and J. Sapirstein, Atomic Data and Nuclear Data 

Tables 64, 279 (1996). 
[32] V.A. Dzuba, V.V. Flambaum, and O.P. Sushkov, Phys. 

Lett. A 142, 373 (1989). 
[33] C.S. Wood, S.C. Bennett, D. Cho, B.P Masterson, J.L. 

Roberts, and C.E. Tanner, Science 275, 1759 (1997). 
[34] S.C. Bennett and C.E. Wieman, Phys. Rev. Lett. 82, 

2484 (1999). 

[35] R. Louden, Quantum Theory of Light, 2nd ed. (Oxford 

University Press, Oxford, 1992). 
[36] C.E. Moore, in Atomic Energy Levels, Natl. Bur. Stand. 

(U.S.) Circ. No. 35 (U.S. GPO, Washington, D.C., 1971). 
[37] Bureau Internationale des Poids et Measures, Pavilion de 

Breteuil, 92312 Svres Cedex, France. 
[38] Th. Udem, J. Reichert, T. Hansen, and M. Kourogi, 

Phys. Rev. A 62, 031801R (2000). 
[39] M.E. Rose, Elementary Theory of Angular Momentum 

(Wiley, New York, 1957). 



